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We study the gauge covariance of the massive fermion propagator in three as well as four dimen- 
sional Quantum Electrodynamics (QED). Starting from its value at the lowest order in perturbation 
theory, we evaluate a non-perturbative expression for it by means of its Landau-Khalatnikov-Fradkin 
(LKF) transformation. We compare the perturbative expansion of our findings with the known one 
loop results and observe perfect agreement upto a gauge parameter independent term, a difference 
permitted by the structure of the LKF transformations. 
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I. INTRODUCTION 



\^ ' In a gauge field theory, Green functions transform in a specific manner under a variation of gauge. In Quantum 
Electrodynamics (QED) these transformations carry the name Landau-Khalatnikov-Fradkin (LKF) transformations, 
. Later, they were derived again by Johnson and Zumino through functional methods, . These transformations 
are non-perturbative in nature and hence have the potential of playing an important role in addressing the problems 
of gauge invariance which plague the strong coupling studies of Schwinger-Dyson equations (SDE). In general, the 
rules governing these transformations are far from simple. The fact that they are written in coordinate space adds 
^v^j . to their complexity. As a result, these transformations have played less significant and practical role in the study of 
SDE than desired. A consequence of gauge covariance are Ward-Green- Takahashi identities (WGTI), |^-^], which are 
simpler to use and, therefore, have been extensively implemented in the above-mentioned studies. 

The LKF transformation for the three-point vertex is complicated and hampers direct extraction of analytical 
restrictions on its structure. Burden and Roberts, |2q| , carried out a numerical analysis to compare the self-consistency 
of various ansatze for the vertex, JlQLp2[, by means of its LKF transformation. In addition to these numerical 
fS |' constraints, indirect analytical insight can be obtained on the non perturbative structure of the vertex by demanding 
correct gauge covariance properties of the fermion propagator. There are numerous works in literature, based upon 
this idea, |L2| - [l7|l Q However, all the work in this direction has been carried out in the context of massless QED3 and 
QED4. The masslessness of the fermions implies that the fermion propagator can be written only in terms of one 
function, the so called wavefunction renormalization, F(p). In order to apply the LKF transform, one needs to know 
. £h ! a Green function at least in one particular gauge. This is a formidable task. However, one can rely on approximations 
based on perturbation theory. It is customary to take F(p) = 1 in the Landau gauge, an approximation justified 
by one loop calculation of the massless fermion propagator in arbitrary dimensions, see for example, |lS[ ]. The LKF 
transformation then implies a power law for F(p) in QED4 and a simple trigonometric function in QED3. To improve 
upon these results, one can take two paths: (i) incorporate the information contained in higher orders of perturbation 
theory and (ii) study the massive theory. As pointed out in jll|, in QED4, the power law structure of the wavefunction 
renormalization remains intact by increasing order of approximation in perturbation theory although the exponent of 
course gets contribution from next to leading logarithms and so on In |l6| , constraint was obtained on the 3-point 
vertex by considering a power law where the exponent of this power law was not restricted only to the one loop 
fermion propagator. In QED3, the two loop fermion propagator was evaluated in [ p|[21p^ ], where it was explicitly 
shown that the the approximation F(p) — 1 is only valid upto one loop, thus violating the transversality condition 
advocated in jl7|]. The result found there was used used in |2^] to find the improved LKF transform. 

In the present article, we calculate the LKF transformed fermion propagator in massive QED3 and QED4. We 
start with the simplest input which corresponds to the lowest order of perturbation theory, i.e., S(p) = m in 

the Landau gauge. On LKF transforming, we find the fermion propagator in an arbitrary covariant gauge. In the 



1 A criticism of the vertex construction in jl^] was raised in jl8) 

2 For the two loop calculation of the fermion propagator, see for example 



1 



case of QED3, we obtain the result in terms of basic functions of momenta. In QED4, the final expression is in the 
form of hypergeometric functions. Coupling a enters as parameter of this transcendental function. A comparison 
with perturbation theory needs the expansion of the hypergeometric function in terms of its parameters. We use the 
technique developed by Moch et. at, |24j], for the said expansion. We compare our results with the one loop expansion 
of the fermion propagator in QED4 and QED3, ]25|,^6|], and find perfect agreement upto terms independent of the 
gauge parameter at one loop, a difference permitted by the structure of the LKF transformations. We believe that 
the incorporation of LKF transformations, along with WGT identities, in the SDE can play a key role in addressing 
the problems of gauge invariance. For example, in the study of the SDE of the fermion propagator, only those 
assumptions should be permissible which keep intact the correct behaviour of the Green functions under the LKF 
transformations, in addition to ensuring that the WGTI is satisfied. It makes it vital to explore how two and three 
point Green functions transform in a gauge covariant fashion. In this article, we consider only a two point function, 
namely, the fermion propagator. 



II. FERMION PROPAGATOR AND THE LKF TRANSFORMATION 



We start by expanding out the fermion propagator, in momentum and coordinate spaces respectively, in its most 
general form as follows : 

S F {x',Z) = fX{ X ',Z)+Y(x;Z), (2) 

where F(p; £) is generally referred to as the wavefunction renormalization and M(p; £) as the mass function. £ is the 
usual covariant gauge parameter. Motivated from the lowest order perturbation theory, we take 

F(p; 0) = 1 and M{p; 0) = m . (3) 

Perturbation theory also reveals that this result continues to hold true to one loop order for the wavefunction renor- 
malization. Eqs. (p]|3) are related to each other through the following Fourier transforms : 

SfM = I d d xe^ x S F (x;0 (4) 



Sf(x;0 = J ^e-^S F (p;0, (5) 

where d is the dimension of space-time. The LKF transformation relating the coordinate space fermion propagator 
in Landau gauge to the one in an arbitrary covariant gauge reads : 

S F {x; £) = S F (x- o) e -^ d (o)-A d (x)] ^ (6) 

where 

Jo (27iT V 

e 2 is the dimensionless electromagnetic coupling. Taking ip to be the angle between x and p, we can write d d p = 
dpp d ~ x sin d ~ 2 ^Q d _ 2 , where il d - 2 = 2 ■ K ( d - 1 )/ 2 /V {{d - l)/2). Hence 

nOO pit 

& d {x) = -ite 2 ^- d f(d) dpp d - 5 / di/jsm d - 2 i;e- ipxco ^ , (8) 
Jo Jo 

where f(d) = '(27r) . Performing angular and radial integrations, we arrive at the following equation 

A d(x) = --^-(^-"r U - 2] . (9) 



With these tools at hand, the procedure now is as follows : 

• Start with the lowest order fermion propagator and Fourier transform it to coordinate space 

• Apply the LKF transformation law. 

• Fourier transform the result back to momentum space. 
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III. THREE DIMENSIONAL CASE 



Employing Eqs. the lowest order three dimensional fermion propagator in Landau gauge in the position 

space is given by 

e- mx (l + mx) , , 

X(x;0) = i-j-^, (10) 



me 

Y{x;0) = 



-nix 



Airx 

Once in the coordinate space, we can apply the LKF transformation law using expression (Jj^) explicitly in three 
dimensions : 

A 3 (x) = -^, (11) 
where a = e 2 /An. The fermion propagator in an arbitrary gauge is then 

S(x;t) = S F (x;0)e-W> x . (12) 
For Fourier transforming back to momentum space, we use 



Performing the angular integration, we get 

nOC 

A(p;t) = / dx smpxe~( m+a ^ x , 

P Jo 

B(p;£) = — f — (1 + mx) [px cos px — sin px] e 
P Jo x 



(14) 
(15) 



and the radial integration then yields 



4m 

Ap 2 + (2m + at) 
„/ ^ Ap 2 + a£(2m + at) at 

B( P ; t) - ~ f 2 1 (o T 7W + 9 arctan 2p / 2m + Q ^ ' 17 
Ap^ + (2m + at) 2p 

One can now arrive at the following expressions for the wavefunction renormalization and the mass function, respec- 
tively : 

at r 2p(4p 2 + a 2 £ 2 )-a£(4p 2 + a£(2TO + a£)) arctan [2p/(2m + a^)] 

(P; ^ ~ ~ 2^ arCtan 1 P/ ( + <)J + 2p(4p 2 + at(2m + at)) - at(Ap 2 + (2m + at) 2 ) arctan [2p/(2m + <)] ' ( } 

£ ) = ^ (19) 

y ' 2p(Ap 2 + at(2m + at)) - at(Ap 2 + (2m + at) 2 ) arctan [2p / (2m + at)} ' 

In the massless limit, one immediately recuperates the well-known results : 

at 2p 

^massless (PS t) = 1 ~ 7T" arctan ~ 7 i ( 20 ) 

2p a£ 

massless (P', t) = . 

In the weak coupling, we can expand out Eqs.( )l1||l9| ) in powers of a. To 0(a), we find 

at 

F(p;t) = 1 H o [(to 2 — p 2 ) arctan [p/ml — mp] , (21) 

2p d 



M (p; t) = to 



1 + — 4r {(to 2 +p 2 ) arctan [p/ml — mp} 
2jr 



(22) 



The expression for the wavefunction renormalization function fully matches with that obtained in pq] , while the one 
for the mass function is also in agreement upto a term proportional to at°, as allowed by the structure of the LKF 
transformations. 
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IV. FOUR DIMENSIONAL CASE 



Employing Eqs. tne l° wes t order four dimensional fermion propagator in position space is given by 

X(x;0) = -— r - 5 K 2 (mx), (23) 

9 

772 

Y{x ] Q) = -— r K^mx) , (24) 

where K\ and K 2 are Bessel functions of the second kind. In order to apply the LKF transformation in four dimensions, 
we expand Eq. (0) around d = 4 — e and use the following identities 



n-2) =-"7 + 0(e), 
x e = l + elnx + C(e 2 ) 

to obtain 

£e 2 



A 4 (x) = i 



167T 2 - e / 2 



- + 7 + 2 In fix + (9(e) 



(25) 



Note that we cannot write a similar expression for A4(0) because of the presence of the term proportional to In a;. 
Therefore, we introduce a cut-off scale x m i n . Now 



where v = at^/Ait. Hence 



x 2 V 

A 4 (x mm ) - A 4 (x) = -iln[ — — j , (26) 

■ rain / 



S F (x;$=S F (x;0)[ — ) . (27) 



For Fourier transforming back to momentum space we use the following expressions 

On carrying out angular integration, we obtain : 

A(p;0 = / dxx- 2 » +1 K^mxJJxipx), (29) 

P Jo 

B(p;0 = -m 2 / dxx~ 2u+1 K 2 (mx) J 2 (px) . (30) 
Jo 

The radial integration then yields : 



A(p;0=~ fe) r(l-i/)r(2-i/) , C!l) 



2 \ v / 1 

m \ A 2 

where we have identified 2/x m i n — ■> A. The above equations imply 
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r(i-iz) 



2m 2 T(3-is) 2 ^i(l-^3-^;3;-p 2 /m 2 ) V A 2 



2\ v 



4m 2 r 2 (2-i/) 2 F 2 M-i/,2-1/; 2;-^ + p'P , (3-i/) 2 F/ 1 - v, 3 - i/; 3; 



2to 2 Fl (l - v, 2 - i/; 2; -p 2 /r 



(2-iz) 2 fi(l-i/,3-i/;3;-p7m 2 ) 



(33) 



(34) 



Eqs. (p3| , p4| ) constitute the LKF transformation of Eqs. (|3j). We shall now see that although Eqs. (|3j) correspond 
to the lowest order propagator, their LKF transformation, Eqs. ( (33"|]34l ) , is non-perturbative in nature and contains 
information of higher orders. 



A. Case a = 



Let us switch off the coupling and put a = which implies v — 0. Now using the identity 
2 F 1 (l,2;2;-p 2 /m 2 ) = 2 F 1 (l,3;3;-p 2 /m 2 ) = (1 + p 2 /m 2 )" 1 , 

it is easy to see that 

F(p;0 = l and M(p;Q=m, 
which coincides with the lowest order perturbative result as expected. 



(35) 
(36) 



B. Case m >> p 

In the limit m >> p, the hypergeometric functions in Eqs. (|33|j34|) can be easily expanded in powers of p 2 /in 2 , using 
the identity 



2-Fi [a, (3; 7; = 1 2 + °\~2 



Retaining only 0{p 2 /m 2 ) terms, we arrive at 



r(i - u)r(2 - v) 



2v ( hv\ p- I p 

3 \ 8 / vn V to 

2 x 2 



2 x -> 

2" 



m 
A 2 " 



1 + 5(1-^ + 

m V m z 



(1-I//2) 

Now carrying out an expansion in a and substituting v — at; /An, we get the following 0(a) expressions : 



F(p;0 = l 



4tt 



27- 
8tt 



1 2p 2 



In- 



2 3m 2 A 2 



1 



P 



3m 2 



(37) 

(38) 
(39) 

(40) 
(41) 



Let us now compare these expressions against the one-loop perturbative evaluation of the massive fermion propagator, 
see e.g., p§ : 



Fi_ l00 p(p;C) = 1 - 

■Mi_i 00 p(f>;£) = m 



ill 

am 

7T 



4/ 4 



£ m 2 



I - 7 ) + 7 (C^-L) + J-^-(l-i) 



(42) 
(43) 
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where 



L = 
C = 



m 2 \ , / p 2 
p z J \ m z 



2 / 771 

7 — In 7r — In — 



Knowing the fermion propagator even in one particular gauge is a prohibitively difficult task. Therefore, Eqs. (j^) 
have to be viewed only as an approximation. For the wavefunction renormalization F(p), this approximation is valid 
upto one loop order, whereas, for the mass function, it is true only to the lowest order. Therefore we cannot expect 
the LKF transform of Eqs. (|3j) to yield correctly each term in the perturbative expansion of the fermion propagator. 
However, it should correctly reproduce all those terms at every order of expansion which vanish in the Landau gauge 
at 0(a) and beyond. Therefore, we expect Eq. (fl2|) to be exactly reproduced and Eq. ( f43|) to be reproduced upto the 
terms which vanish in the Landau gauge at 0(a). After subtracting these terms, we call the resulting function as the 
subtracted mass function : 



loop 



(p;0 



a£m 
Air 



1 



4o- 



L) 



In the limit 



oo, the wavefunction renormalization acquires the form 



a£ 

F(p;Oi-loop = l + ^ 



2p 2 

377J 2 



while the subtracted mass function is 



at, 
8^ 



3m 2 



(44) 



(45) 



(46) 



The last two expressions are in perfect agreement with Eqs. (^,^l|) after we make the identification 



C/i e -> 2 7 + ln 



A 2 



(47) 



C. Case of weak coupling 



The case m >> p is relatively easier to handle as we merely have to expand 2-Fi(/3, 7; 5; x) in powers of x and retain 
only the leading terms. If we want to obtain a series in powers of the coupling alone, we need the expansion of the 
hypergeometric functions in terms of its parameters (3 and 7. We follow the technique developed in p4| . One of the 
mathematical objects we shall use for such an expansion are the Z-snvas defined as : 



Z(n;mi, . . . ,mk;xi, ...,Xk) 



ra>ii>i 2 >->ifc>0 



mi 
1 



X fc 
■nifc 



For xi = . . . = Xk = 1 the definition reduces to the Euler-Zagier sums, 27 2q] : 

Z(n;mi,...,m k ;l, ■■■,!) = Z mit ... :mk {n) . 
Euler-Zagier sums can be used in the expansion of Gamma functions. For positive integers n we have 

r(n + e) = r(l + e)r(n) [l + eZ x {n - I) + . . . + e n-1 Z u .„ 1 (n - 1)] . 
The first sum Z\[n — I), e.g., is just the (n — l)-th harmonic number, i?„_i, of order 1 : 



n — 1 

Z 1 (n-l) = y2- = H n _ 1 



(48) 

(49) 
(50) 

(51) 



G 



With these definitions in hand, we proceed to expand a hypergeometric function, 2-Fi(l +e,2 + e;2;x), as an example, 
assuming |x| < 1 : 



2 F 1 {l + e,2 + e:2;x) = 1 



= 1 



r(2) 



n i. i :)rC2 he) 

v ' • J n— 1 

1 

(l + e)r 2 (l + e) 



r(l + e + n)r(2 + £ + n) a:" 
r(2 + n) n\ 

(1 + £ + n)(e + n) 2 T 2 (e + n) a;" 

r(2 + n) ~n\ 



Employing Eq. ( |50| ) , we can expand the last expression in powers of e to any desired order of approximation. We shall 
be interested only in terms upto 0{a). 



2 F 1 (l+e,2 + e:2;x) = 1 + ]T x n - e ]T x» + £ £ 



2 + 3n 
J n(n + 1) 



2£ ^ F„_ 



(52) 



Performing the summations, we obtain 



2 i ;i i(l + £,2 + £;2;a;) = 



l-x 



■ 1 + X , . 

1-£<M + ln(l-a; 

1 a; 



(53) 



Similarly, 



2^1(1 + 6,3 + 6; 3; a:) 



1 - x 



1 3 1 



a; 2 1 — a; 



1 + x 



l-x 



In (l-x) 



Substituting back into Eqs. (|33|j34| ) and identifying e = —u, we obtain 



(54) 



_2 7 -l n _+ M_ (1-1 



M(p;£) = m + 



47T 



A 

p 2 



(55) 



which matches exactly onto the one loop result of Eqs. ( |42| . |44D after the same identification as before, i.e., (|47 
Therefore, we have seen that the LKF transformation of the bare propagator contains important information of 
higher orders in perturbation theory. 



V. CONCLUSIONS 



We have studied the gauge covariance of the massive fermion propagator in three as well as four dimensional 
QED through its LKF transformation, starting from its lowest order approximation. Eqs. (|l^,|l^,^,Q) form the 
main result of this article. In the three dimensional case, the LKF transformation consists of basic functions of 
the momentum variable, whereas, in the four dimensional case, hypergeometric functions arise with electromagnetic 
coupling as parameter of these functions. Although our input is only the bare propagator, the corresponding LKF 
transformation, being non-perturbative in nature, contains useful information of higher orders in perturbation theory. 
For example, we have shown that a perturbative expansion of our results matches onto the known 1-loop results 
upto gauge independent terms at this order. This slight difference arises due to our approximated input and can be 
corrected systematically at the cost of increasing complexity of the integrals involved. We intend to carry out similar 
exercise for the 3-point fermion-boson vertex. LKF transformations of the propagator and the vertex impose useful 
constraints on the SDE and we believe that these transformations can be of immense help in addressing the problems 
of gauge invariance in the related studies. 
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APPENDIX 



Most of the integrals involved in this paper are listed below for a quick reference 

fd-l 



J dip sin d ~ 2 i/; cos ^e'^ cos * = -i ^ ( y J 2 T 

x J d/2 (ax) = 2l _ d/2 ad/2 
For the three dimensional case, the needed integrals are : 

dO sin 9 e- ipxcos6 



Jd (px) , 



2 sin pa; 



px 



dp 



d6 cos9 sin (9 e - ipxcose = 2i 
cos px sin px 



o (p 2 +m 2 ) 



px (px) 2 



cospx smpx 
px (px) 2 

tt (1 + mx) ™ ; 



dp 



psmpx 7T 



' (p 2 + m 2 ) 2 

1 f°° dx a p 

— / — TT-e [px cospx — smpxj = —1 H — arctan— , 

PJo ■' ' P " 

1 f°° dx 

— / — e ax \px cospx — sinpxl 
P Jo x 

dxsmpxe- {m+alL/2)x 



a 2 + p 2 p 



i arctan^ 
a 



P 

(m + a£/2) 2 +p 2 

For the four dimensional case, we used the following integrals in particular : 



desm 2 ee- ipxcost 

p» +1 J u (px) 



px 



Ji(px) , 



dp 



-Kv-^mx) 



o ^(p 2 + m 2 )^+i 2/T(/i + l)"" p 

A -"- 1 fe ,y /i/-A + /x + l\ (v-X-pi + l 



dxx K^(ax)J v (bx) 



2 x + 1 T(l + v) 

x 2^1 



— A + /i + 1 ^ — A — /i + 1 



Some of the series used in our calculation are as follows 



H n-lX r ' 

n=l 

V H+1 x r. 

^ n(n + 2) 

n— 1 v ' 



x ln(l — x) 
1-x ' 

2 + x (l + x 2 )ln(l-x) 



4x 



2x 2 
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